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Abstract: Recently the ATLAS collaboration reported a Sa excess in the leptonic- 
Z-\-jets-\-E'!p'‘^^ channel. This may be interpreted in the Next-to-Minimal Supersymmetric 
Standard Model (NMSSM) by gluino pair production with the decay chain g —)• qqx 2 
qqZx^, where X 2 denote the lightest and the next-to-lightest neutralinos with 

singlino and bino as their dominant components respectively. After exploring the relevant 
parameter space of the NMSSM by considering the constraints from the ATLAS searches 
for jets + signals, we conclude that the NMSSM is able to explain the excess at la 

level with the number of the signal events reaching its measured central value in optimal 
cases, and the best explanation comes from a compressed spectrum such as my ~ 650GeV, 
m^o ~ 565GeV and m^o ~ 465GeV. We also check the consistency of the ATLAS results 
with the null result of the CMS on-Z search. We find that under the CMS limits at 95% 
C.L., the event number of the ATLAS on-Z signal can still reach 11 in our scenario, which 
is about 1.2a away from the measured central value. 
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1 Introduction 


Since the discovery of a Higgs-like particle in 2012 [1], the search for new physics beyond 
the Standard Model (SM) has been the most important task of the LHC. So far this 
search has covered a wide range of possible signatures of new physics, especially the signals 
of supersymmetry (SUSY), which are characterized by various combinations of missing 
transverse energy plus jets and/or leptons (electrons or muons). In this direction, 

one impressive result is that the ATLAS collaboration recently reported a 3a excess in the 
channel of two leptons with an invariant mass located around mz, large and at least 

two jets [2]. For 20.3/6“^ of integrated luminosity at a center-of-mass energy of 8TeV, the 
number of the observed leptonic-Z events is 29 (summing over electron and muon pairs) 
in comparison with 10.6 ± 3.2 events expected in the SM [2]. 

During past weeks, several works appeared to interpret the excess in SUSY [3-7] and all 
of them employed gluino pair production process. As pointed out in [3], in order to enhance 
the signal rate to meet the ATLAS data, the gluino must be relatively light, nig < 1.2TeV, 
and meanwhile have a rich Z-boson yield in its decay. The latter requirement can be 
satisfied in the General Gauge Mediation model (GGM) with a very light gravitino G as 
the lightest supersymmetric particle (LSP) [8]. In this framework, due to the very weak 
couplings of G to matter fields, gluino prefers to decay firstly into a neutralino, and the 
subsequent decay of the neutralino into G is able to produce one or more Z-bosons [3]. 
The simplest decay chain can be written as 5 —>■ qqxi qqZG, where x\ denotes the 
lightest neutralino and the first step decay proceeds via exchanging an off-shell squark. It 
should be noted that, since nig 3> mz, one has following mass relations: mg ^ m^o > mz 
or mg > m^o 3> mz- In either case, the gluino cascade decay will produce at least two 
hard jets, which are from the first step decay g qqXi or from the hadronic decay of an 
energetic Z boson. As a result, this explanation could be tightly constrained by the search 
for SUSY in events with jets plus Moreover, as illustrated in [4], this explanation can 

not reproduce the shape of the distribution measured by the ATLAS collaboration. 

Given the deficiencies of the GGM in explaining the excess, we in this work con¬ 
sider an alternative simple scenario realized in the Next-to Minimal Supersymmetric Stan¬ 
dard Model (NMSSM) [10] . Our scenario involves a singlino-dominated LSP Xi a-iid a 
bino-dominated next-to-lightest supersymmetric particle (NLSP) X 21 decay chain 

g —7- qqx 2 QQ^Xi- Similar to the GGM explanation, the LSP in our scenario has very 
weak couplings to matter fields, but different from a very light G in GGM (which is re¬ 
quired by the prompt decay of the NLSP [9]), now the LSP may be massive. This feature 
enables us to choose a slightly compressed SUSY spectrum mg ~ m^o ~ m^o + mz to 
evade the constraints from the SUSY searches. Gonsequently, gluino as light as TOOGeV is 
still experimentally allowed, which is helpful to explain the excess (see discussion below). 
Another advantage of our scenario is that it can reproduce the and Ht distributions 

measured by the ATLAS collaboration. This has been recently emphasized in [4]. 

About our scenario, two points should be noted. One is that, since Z boson de¬ 
cays dominantly hadronically, in most cases the gluino cascade will produce multiple jets. 
Gonsidering that the ATLAS collaboration has performed several comprehensive searches 
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for the jets + E'™*® signals, and each of the searches focused on multiple signal regions 
which are sensitive to different configurations of SUSY spectrum [17-19], one should con¬ 
sider these searches by reproducing them through detailed simulations and then using the 
corresponding experimental data as input to limit SUSY. Studying the capability of our 
scenario under these constraints to interpret the excess is the main aim of this work. The 
other is that the CMS collaboration recently also reported its analysis on the leptonic- 
Z -|- jets + signal [20], which should have the same physical origin as the ATLAS 

excess in our scenario. Since the CMS observed no excess, one should check the consistency 
of the two searches by noting that they are based on different cuts and different detectors. 
This is another aim of this work. 

This work is organized as follows. In Section II, we introduce our scenario and illustrate 
its key features. Then in Section III, we describe the strategies of the ATLAS and CMS 
searches for the leptonic-Z + jets + signal as well as other searches for jets + E™*^^ 

signals. In Section IV, we first show the dependence of various signals on SUSY parameters 
in our scenario. Then we scan the relevant SUSY parameters by considering the constraints 
from the SUSY searches, and in the surviving parameter space, we investigate the capability 
of our scenario in explaining the ATLAS Z-peaked excess and check the consistency of the 
measurements from the two collaborations on the leptonic-Z-|-jets + E™'^^ signal. Finally, 
we draw our conclusion in Section V. The validations of our simulations on the ATLAS 
searches for the leptonic-Z + jets + E™*^ signal are presented in the appendix. 

2 A potential explanation of the excess in the NMSSM 

The NMSSM is the simplest singlet extension of the Minimal Supersymmetric Standard 
Model (MSSM) and its superpotential is given by [10] 

Wnmssm = Wmssm + • Hd + CfS + ( 2 . 1 ) 

where Wmssm is the superpotential of the MSSM, and Hd are MSSM Higgs super¬ 
fields, 5 is a gauge singlet superfield, A and k are dimensionless Yukawa couplings, fi' is 
a supersymmetric mass and with dimension of squared mass parameterizes a tadpole 
term. 

Due to the addition of the singlet field, the NMSSM predicts three CP-even Higgs 
bosons, two CP-odd Higgs bosons and five neutralinos. The mass matrix of the neutralinos 
in the bases —iW^, H^, H^, S^) is given by[10] 
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where = cos 9w, Mi and M 2 are soft gaugino masses, Vu = usin/3 and Vd = vcosfj are 
vacuum expectation values (vev) of the fields Hu and Hd respectively, and /Xe// = fJ- + Xs 
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with s denoting the vev of the singlet scalar field S. This matrix indicates that the mixings 
between the bases depend on the parameters tan and A, and for a singlino-dominated 
LSP, its gaugino components should be very small, while its Higgsino components may be 
sizeable for a large A. 

The interactions of the neutralinos take following form [10] 
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where N is the rotation matrix to diagonalize the neutralino mass matrix with its element 
denoted by Aijj {i,j = !,••• ,5), 0!^j = -Ofj* = -^NisN*^ + lNi 4 N*^, and Ilfj' = NiaNjb + 
dXibdXja- This Lagrangian indicates that the interactions of the neutralinos with a light 
quark and those with a Z-boson are determined by their gaugino components and Higgsino 
components respectively, while their interactions with a Higgs boson are decided by the 
both and there may exist cancelation between the two contributions. These features are 
helpful in understanding our scenario. 

As mentioned in last section, the key issues in explaining the leptonic-Z excess are 
how to improve the Z yield rate in gluino decay and how to escape the constraints from 
the ATLAS direct searches for SUSY. Our scheme is as follows: 


• We set the LSP to be the singlino-dominated neutralino, the NLSP to be the bino- 

dominated neutralino, and M 2 ,yeff > iRg- This setup guarantees that the rate of 
the gluino decay into the X 2 intermediate off-shell squark) is much larger than 

that of the decay into the Xv Meanwhile, since the NLSP is not accompanied by a 
light chargino, the gluino may decay into the X 2 ^ approximate 100%. 

• We assume m^, S> rrig > rrig where q denotes light quarks. This assumption 

ensures that the gluino decays dominantly into qqx 2 state at hrst step. We suppress 
the decay chain g —)• tix 2 tt^Xi due to its very low efficiency in getting the desired 

signal (We obtain this conclusion by detailed simulations). We assume a relatively 
heavy q since it is favored by current LHC search for SUSY. 

As will be shown later, we are particularly interested in the situation nig ~ In 
this case, the assumption on the mass spectrum may be relaxed to be m^, nig > nig 
and nig < ni^o + 2mt. 

• We require that the decay X 2 ^Xi proceeds at a rate of approximate 100%. 
This can be realized by setting mz < m^o — m^o < nih or by choosing specific 
SUSY parameters such that the coupling hxiX 2 suppressed. With the package 
NMSSMTools[ll], we numerically checked that, by tuning the parameters of the 
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NMSSM, the rate of X 2 ^Xi can be significantly larger than that of X 2 ^Xi 
even when m^o — m^o > mh- 

• We employ a moderately compressed spectrum nig ~ m^o ~ m^o + mz to avoid the 
constraints from the LHC searches for SUSY. 

About our scenario, we have more explanations. One is that the measured SM-like 
Higgs boson mass in our scenario can be satisfied by choosing appropriate values of A, tan (3 
and also third generation squark masses [12-15] and the correct dark matter relic density 
can be achieved by the annihilate channels XiA? Hi/Aj —)• XY with XY being any 
possible final states [16]. The other is that, in our scenario, we ad hoc forbid the possibility 
of M2, < ■'TT-g- The main reason is that in that case, the gluino decay chain may become 

more complicated, and channels such as g ^ Qq'xf are open. 

Then one has to scrutinize the vast parameter space of the NMSSM to find parameter 
points where the Z yield rate is high. Since this process involves the simulation of the 
gluino pair production at the LHC, it is very time-consuming in calculation by cluster 
because the decay chains are lengthy. So our scenario should be regarded as a simplest 
attempt to explain the excess, which only involves the parameters m^o and m^o. 

3 Simulations and cuts 

In order to discuss the constraints of the LHC searches for SUSY on our scenario and 
also calculate the event number of the leptonic-Z signal, we implement the corre¬ 

sponding experimental searches in the package CheckMATE-1.2.0 [23], which have tuned 
the package DelphesS.O.lO [24] to simulate the behavior of the ATLAS (CMS) detector. 
The validations of our simulations are presented in Appendix for the benchmark points 
provided by the experimental groups, and we find that our results agree with the corre¬ 
sponding experimental analyses at 20% level. Throughout this work, the events in our 
simulation are generated by MG5_aMC[25] which includes Pythia[26] for parton shower¬ 
ing and hadronization, and the package Prospino [27] is used to calculate the NLO cross 
sections of the gluino pair production. 

In the following, we only briefly describe the cut flows of these experiments. Detailed 
information about them can be found from the corresponding experimental reports. 

3.1 ATLAS search for leptonic-Z + jets + signal 

This search concentrated on the events with a same flavor opposite-sign (SFOS) dilepton 
pair, jets and at 8-TeV LHC [2]. One intriguing result of this search comes from 

the signal region SR-Z where the invariant mass of the lepton pair locates around mz- 
Explicitly speaking, in contrast with 4.2 ± 1.6 electron pair events and 6.4 ± 2.2 muon pair 
events expected for the SM background, 16 events and 13 events are observed respectively 
in this region. This implies a 3cr excess of the observed events over the background, and 
may be regarded as a hint of SUSY. So in the following we try to explain this excess with 
our scenario. The signal region SR-Z is defined by 
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Table 1: Definition of signal regions in ATLAS search on > 6 jets + signal [19]. 

Here the jet, the pseudo-rapidity r] and multiplicity all refer to the R = 0.4 jets. Composite 
jets with a larger radius parameter i? = 1.0 are used in the multi-jet -|- Mj stream when 
constructing the quantity Mj (see text). 



Multi-jet + flavour stream 

Multi-jet J- m'^ stream 

Name 

8j50 

9j50 

>10j50 

7j80 

>8j80 

>8j50 

>9j50 

>10j50 

Jet P™"[GeV] 

50 

80 

50 

Jet count 

=8 

=9 

>10 

=7 

>8 

>8 

>9 

>10 

b-jet count 

0,1,>2 

0,1,>2 


0,1,>2 

0,1,>2 


M^lGeY] 


> 340, > 420 


• The first two leading leptons are SFOS, their transverse momentums are larger 
than 25GeV and lOGeV respectively, and meanwhile their invariant mass satisfies 
SlGeV < mu < lOlGeV. 

• More than one jets have Pt > 35GeH and \r]\ < 2.5. 

• > 22bGeV and Ht > GOOGeV, where Ht is defined as the scalar sum of the 
transverse momenta for all signal jets and the two leading leptons. 

• The azimuthal opening angle between each of the leading two jets and should 

be larger than 0.4, i.e. A4>{jeti^2, > 0.4. 

3.2 ATLAS search for > 6 jets + signal 

To search for final states with large jet multiplicities and ATLAS used two inde- 

y^ 

pendent streams, i.e. multi-jet -|- flavor stream and multi-jets-|-Mj stream, in its anlaysis 
[19]. 

In the multi-jet -|- flavour stream, the number of jets with \r]\ < 2 and px above a 
threshold (50GeV or SOGeV) is firstly determined for a given event. Then events are 
categorized by the number of jets, and further subdivided by the number of jets that are 
6 -tagged. After such treatments, totally thirteen signal regions are defined, which are given 
in Table. 1. 

y^ 

Analysis of the multi-jet -|- Mj stream seems more complicated. As the first step, the 
number of {R = 0.4) jets with |r/| < 2.8 and px above 50GeV is determined for a given 
event, and like the former stream, events are classified by their jet numbers. As the second 
step, the four-momenta of the R = 0.4 jets are used as inputs for a second iteration of 
the anti-fc* jet algorithm [21], which now adopts a larger jet radius parameter, R = 1.0. 

y^ 

The resulting larger objects are called as composite jets. Then a selection variable Mj is 
defined by 

= (3.1) 

j 
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where the sum is over the mass of the composite jets, and two signal regions are defined 
by the threshold of Mj (340GeV or 420GeV). In this stream, there are totally six signal 
regions (see Table. 1). 

We mention by the way that both the streams veto events containing any isolated 
electron or muon candidates with Pt > IQGeV at the start of the analysis, and require 
EJp^^^/y/IPp > 4GeI4i/2 at the end of the cut flows, where Ht is the scalar sum of the 
transverse momenta of all jets with Pt > 4:0GeV. 

3.3 ATLAS search for 2 ~ 6 jets + signal 

The search for events with 2 ~ 6 jets + Etp^^^ at 8 -TeV LHG with a total integrated 
luminosity of 20.3/6“^ has been carried out in [17] and also in [18]. Gompared with the 
preliminary search described in [17], the updated search in [18] defined more signal regions 
to optimize the search. For example, it has been shown that the signal regions 2jW and 
4:jW are able to improve greatly the sensitivity of the ATLAS search to models predicting 
the cascade g —)• qq'xf Q^^Xi- The main cuts of this search are: 

• > leOGeV. 

• -Pr(ji) > 130Gel/ for the leading jet, and PtU) > OOGeV for the other jets. An 
exception is the signal region 4jlT where PtU) > ^OGeV is required. 

• The azimuthal angles between jets and the satisfy 64>{jeti 2 ,{ 3 )y > 0.4 

and 64 >{jeti^ 3 ,Etp^^^)min > 0 . 2 . 

• A variety of signal regions is then defined by the number of jets, the values of 

Ejp^^^/^H t-, Etp^^^/nieffiNj) and meff{incl). Here Ht, meff{Nj) and meff{incl.) 
are defined as the scalar sum of the transverse momenta for all Pt > AOGeV jets, the 
leading Nj jets and and all jets with Pt > 40GeV and EJp^^^ respectively. 

In practice, we also consider the search in [17]. This search has been implemented in 
the GheckMATE by the package authors, and we have verified its correctness. Since this 
search is significantly weaker than that of [18] in constraining SUSY parameter space, we 
do not list its cuts here. 

3.4 CMS search for leptonic-Z Pjets + Etp^^^ signal 

In [20], the GMS collaboration carried out a dedicated on-Z counting experiment to search 
for leptonic-Z +jets + Fi™*^ signal. This signal is same as that of the ATLAS search in 
[ 2 ], but different from the ATLAS result, the GMS collaboration saw no excess. Since the 
two signals have same physical origin in our scenario, they should be correlated. So in our 
work, we consider the CMS search as a possible constraint on our interpretation of the 
excess. The main strategies of the CMS search are given by 

• There exists at least one SFOS dilepton pair (electron or muon pair) with Pt > 
20GeV and \r]\ < 2.4 for each lepton, and the invariant mass of the lepton pair 
satisfies 81 < mu < lOlGeU. 
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Table 2: Six bins of the CMS dedicated counting experiment for events with an on-shell 
Z boson [20]. For each bin, the 95% CL upper limits on the number of signal events, 
are also presented. 


^jets 

> 2 

> 3 

Ef*®®[GeU] 

100-200 

200-300 

>300 

100-200 

200-300 

>300 

^data 

1187 

65 

7 

490 

35 

6 

^bkg 

1204±106 

74.5±11.3 

12.8±4.3 

478±43 

39.2±6.6 

5.3±2.3 

c95 

‘^obs 

207 

20 

7.6 

89 

16.1 

8 


• Six bins are defined by the value of and the number of jets with Pt > 40GeV 

and \r]\ < 3.0. Here we present some information of these bins in TABLE 2. 

About the CMS search, we note that, unlike above ATLAS searches, it does not provide 
the 95% C.L. upper limits on the number of signal events for each bin, which was denoted 
as in literature [23]. In order to discuss the constraint of such a search, we calculate 
these limits by the asymptotic CLg prescription [ 22 ], and present the corresponding results 
in TABLE 2. As we will show below, the bins with Nj >2,3 and 200 < < 300 

usually put the strongest limitation on any explanation of the excess, which is different 
from previous analysis [4]. 

4 Numerical results and discussions 

As we introduced in section II, the branching ratio of 5 —)• qqx 2 QQ^Xi our scenario can 
approach 100% by setting the SUSY parameters other than the masses rrig, m^o and m^o at 
certain values. In this case, the event number of the gluino pair production with a certain 
final state is solely determined by the masses, or equivalently by rrig, Ami = 
and Am 2 = m^o — m^o. In this section, we investigate the constraints of the SUSY 
searches on the masses. For this end, we define for each experimental search the ratio 
R = max{Ns^i /where Ns^i is the signal event number for i signal region in the 
search, is its 95% upper limit, and the max is over all signal regions of the search. 

Obviously, only in the case of i? < 1, the corresponding mass spectrum is allowed by the 
search at 95% C.L.. We also investigate how well our scenario can explain the leptonic-E 
excess. After considering the statistical uncertainty of the measured event number, we infer 
from the ATLAS search that the needed signal event number, Nu, is 18.4 ± 6.3 (assuming 
the experimental data obey Poisson Distribution and adding the statistical uncertainty and 
systemic uncertainty in quadrature), which corresponds to 12.1 < Nu < 24.7 and 5.8 < 
A^ii < 31 in explaining the excess at la and 2 ct levels respectively. In following discussion, 
it is convenient to define the normalized leptonic-Z event number as Ru = Nu/18.4, then 
the Icr and 2 ct ranges are restated as 0.66 < Ru < 1.34 and 0.32 < Ru < 1.68 respectively. 

In Fig.l, we show the dependence of Rs for different searches on Ami, Am 2 and mg 
respectively by fixing mq = 4.5 TeV. From this figure, one can learn at least five facts. 
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Figure 1: The dependence of Rs for different searches on Ami, Am 2 and mg. Here 
Ami = mg — m^o denotes the mass splitting between g and X 2 Am 2 = m^o — m^o is 
that between X 2 Xi- For the ATLAS leptonic-Z analysis, R is the normalized lepton 
pair event number, R = Nu/IS.A with 18.4 representing the central value of the excess, 
while for any of the other searches, R = max{Ns^i /where Ns^i is the event number 
for i signal region in the search, is its 95% upper limit, and the max is over all signal 
regions in the search. In the latter case, R < 1 means that the parameter point is allowed 
by the corresponding search at 95% C.L.. 

1. The left and middle panels indicate that, as the SUSY spectrum becomes compressed, 
the value of R for each search decreases monotonously. This reflects that a smaller 
mass splitting is able to suppress the cut efficiency in the search so that fewer SUSY 
signal events are finally retained. Consequently, this case is apt to escape the con¬ 
straints from the direct search for SUSY at the LHC, but meanwhile, the ATLAS 
leptonic-Z events are suppressed too. 

2. Motivated by above observation, we consider a rather compressed spectrum in the 
right panel, and study the dependence of R on mg. The panel indicates that Rs 
decrease monotonously as the gluino becomes heavy. This is because the suppression 
of the gluino pair production rate with the increase of mg usually dominates over the 
enhancement of the cut efficiency in getting the signal events of the search, and the 
net result is the dropping of the signal events. 

3. All the three panels indicate that, in the case of a very compressed spectrum, the 

CMS search for the leptonic-Z +jets -|- signal puts the severest constraint on 

SUSY among the searches. Especially, as the gluino becomes light, the signal event 
number in the CMS search grew much faster than those of the other searches. As a 
result, the CMS search may be used to set a lower bound on mg in our scenario. 
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4. In general, the ATLAS preliminary search for 2 ~ 6 jets + signal and that for 

> 6 jets + signal are comparable in constraining the parameter space of our 

scenario, and the both constraints are weaker than that from the ATLAS updated 
search for 2 ~ 6 jets + signal. 

5. Since the values of Rs for the experimental searches have different behaviors with 
the variation of the mass spectrum, it is better to scan the masses in studying the 
capability of our scenario to explain the excess. 

In Fig. 2, we show the constant contours of the event number for the ATLAS search for 
the leptonic-Z +jets + signal on the Ami — Am 2 plane for different gluino masses. 

Regions between line 12.1 and line 24.7 and those between line 5.8 and line 31 are able 
to explain the excess at Icr and 2a levels respectively. Constraints from the other SUSY 
searches are also plotted by different types of lines with the left side of the lines being 
experimentally allowed. Solid line, dashed line and dotted line represent the constraints 
from the ATLAS searches for > 6 jets + 2 ~ 6 jets + (updated) and 2 ~ 6 

jets + (preliminary) respectively, and they are obtained by fixing the corresponding 

R value at 1. The CMS constraint is plotted in a similar way, but by dash-dotted lines. In 
getting this figure, we scan the parameters Ami and Am 2 with a grid of 30GeV for Ami 
and lOGeV for Am 2 . 

From Fig. 2, one can learn following facts: 

• With the decrease of the gluino mass, only a more compressed spectrum can survive 
the constraints from the LHC searches for SUSY. Among the constraints, the CMS 
search for the letponic-Z +jets + EJp^^^ becomes the most stringent constraint for 
rrig < 800GeV. Especially, in the case of rrig = 700GeV, only a very small corner on 
the Ami ~ Am 2 plane survives the constraint. 

• If one does not consider the CMS constraint, our scenario can explain well the ATLAS 
leptonic-Z excess without conflicting with the ATLAS searches for jets + Etp^^^. The 
central value of the excess (18.4 events) may be achieved for rrig < 700GeV, and one 
of the best benchmark points is rrig = 650GeV, m^o = 565GeV and m^o = 465GeV. 
While if one take the CMS constraint serious, we find that the leptonic-Z events in 
the ATLAS experiment is at most 11. This tension arises since the CMS signal and 
the ATLAS excess have same physical origin. 
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Figure 2: Constant contours of the event number for the ATLAS search for the leptonic- 
Z +jets + signal on the Ami ~ Am 2 plane for different gluino masses. Regions 

between line 12.1 and line 24.7 and those between line 5.8 and line 31 are able to explain 
the excess at Icr and 2a levels respectively. Constraints from the other SUSY searches are 
also plotted by different types of lines with the left side of the lines being experimental 
allowed. For the case of nig < 650GeV, there are no regions surviving the constraint from 
the CMS search for the leptonic-Z signal. 
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5 Conclusion 


Recently the ATLAS collaboration reported a 3cr excess in the leptonic-Z + jets + 
channel. In this work, we intend to explain the excess by the gluino pair production 
process in the NMSSM. For this purpose, we implemented this experimental analysis (and 
also some other useful experimental analyses) in the package CheckMATE, and verihed the 
reliability of our simulations. 

The scenario we considered involves a singlino-like LSP, a bino-like NLSP and an 
approximate 100% branching ratio for the decay chain g —)• qqx 2 QQXi- this case, 
the SUSY signal rates after cuts are essentially determined by the masses rrig, m^o and 
m^o. We scanned these masses by considering the constraints from the ATLAS searches 
for jets + signal, then we concentrated on the leptonic-Z signal in the experimentally 

allowed parameter space. We concluded that our scenario is able to explain the excess at 
la level with the number of the signal events reaching its experimental measured central 
value in optimal cases, and the best explanation comes from a compressed spectrum such 
as Trig ~ 650GeV, m^o ~ 565GeV and m^o ~ 465GeV. We also tested the consistency of 
the ATLAS result with the CMS result in the same channel, which observed no excess on 
the leptonic-Z production. We found that under the CMS limits at 95% C.L., the number 
of the ATLAS signal events can still reach 11 in our scenario, which is about 1.2a away 
from the measured central value. 

We emphasize that, although we employ a specihc scenario in our work to explain 
the ATLAS excess, our conclusions seem to be general. That is, the NMSSM is able to 
explain quite well the excess if only the constraints from the ATLAS searches for SUSY 
are considered. While on the other hand, if one considers the CMS search limits, there is 
a mild tension between the ATLAS excess and the CMS result. 

Note added: At the final stage of this work, there appeared a paper by Ulrich 
Ellwanger [4], which considered a similar scenario to ours. Anyhow, our work differs from 
his work mainly in following three aspects: 

• In his paper, only two benchmark points were discussed, while we perform a scan over 
the relevant parameter space. Our calculation by cluster is very time-consuming. 

• We considered more constraints than in his paper, where only the ATLAS preliminary 

search for 2 ~ 6 jets + was considered. Especially, the benchmark point P2 

in his paper was considered as an optimal point to explain the excess, while in our 
work, we got R = 1.19 for the updated ATLAS search for 2 ~ 6 jets + (with 

rriq = 4:.5TeV). Considering that we have about 20% uncertainty in simulating 
thesearch, we infer that the point is at the boundary of being excluded. 

• When we consider the CMS analysis, we found that the most stringent constraint 

comes from the bins of Njets ^2,3 and 200GeV < < 300GeV for most cases, 

instead of the bin Njets ^ 2 and E™®® > 300GeV adopted in [4]. This is because 
we are considering a compressed spectrum, where the E™®® of the signal is usually 
small. 
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A Validation of Monte Carlo simulations 

In this appendix, we show the validations of our simulations on the ATLAS searches for 
the leptonic-Z +jets + signal [2] and the 2 ~ 6 jets + signal [18]. For each 

search, we choose two benchmark points adopted by the ATLAS report, and repeat the 
experimental analysis. Our simulation results are presented in Table 3 and Table 4 respec¬ 
tively, and all of them indicate that our simulation agrees with the corresponding ATLAS 
analysis within 20% uncertainty. We also validated our simulation on the CMS search for 
the leptonic-Z -\-jets + signal [20] with the corresponding results given in Table 5. 

This table indicates that we can reproduce the CMS results within 10% uncertainty. 

Since the ATLAS preliminary search for 2 ~ 6 jets -|- signal [17] and the search 

for the > 6 jets -|- F™*^ signal [19] do not affect our conclusions, we do not present their 
validation here. In fact, we checked that we can reproduce these searches within 15% 
uncertainties. 
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Table 4: Validation table of our implement of the ATLAS search for the 2 ~ 6 jets + 
CheckMATE. We considered the process pp —)■ gg ^ xfxfjjjj 
XiXi^^W^jjjj at 8-TeV LHC in the MSSM with two benchmark points taken from [18]. 
Numbers in the columns EXP and OUR are the cut efficiencies in 2jW or 6 jt signal region 
obtained by the ATLAS collaboration and us respectively, and those in the column DIEE 
represent their relative difference. 


g g one step, rrig = 1200GeU, m-± = 1150GeU, m^o = GOGeV 

SR;2jW 

EXP 

OUR 

DIEE 

Ejuiss > iQOGeV, Pt(ji,2) > 130(60)GeU 

52.70 

56.35 

-7% 

>0.4 

46.30 

49.01 

- 6 % 

N(W) unresolved > 2 

9.20 

8.70 

5% 

Er^7me//(iV,) > 0.25 

7.00 

6.69 

4% 

meff{incl.) > 1800GeU 

5.30 

4.86 

8 % 

g g one step, rrig = 1265GeU, m-± = 945GeU, m^o = 625GeU 

SR: 6 jt 

EXP 

OUR 

DIEE 

> 160GeU, RrO'i, 7) > 130(60)GeU 

53.30 

54.16 

- 2 % 

PtUs) > 60GeU 

53.00 

53.83 

- 2 % 

Pt^Ja) > 60GeU 

50.50 

51.50 

- 2 % 

PrUb) > 60GeU 

41.40 

43.49 

-5% 

PTiJe) > 60GeU 

26.70 

29.78 

- 12 % 

A7ji,2,3,^r") >0.4 

22.40 

25.38 

-13% 

> 0.2 

18.20 

20.54 

-13% 

EJp^^^/meff{Nj) > 0.25 

10.90 

11.64 

-7% 

meff{incl.) > 1500GeU 

4.20 

4.81 

-15% 
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Table 5: Validation table of our implement of the CMS search for the leptonic-Z +jets + 
jj^miss gigj^al [20] in CheckMATE. We considered the process pp ^ gg ^ XiXiJjjj 
GGZZjjjj at 8 -TeV LHC in the GMSB model with the benchmark point taken from [20]. 
Numbers in the columns EXP and OUR are the event numbers of different bins obtained 
by the CMS collaboration and us respectively, and those in the column DIFF represent 
their relative difference. 



EXP 

OUR 

DIFF 

All events 

37.7 

37.7 

0 

> 2 leptons(/^/'^),PT > 20 GeU 

11.9 

11.8 

-0.5% 

81 < mu < lOlGeU 

10.7 

10.4 

-2.7% 


^jets ^ 2 



EXP 

OUR 

DIFF 

EXP 

OUR 

DIFF 

‘^jets ^ 2 OT 3 

10.7 

10.4 

-2.9% 

10.4 

10.2 

-1.9% 

Ejniss > looGeU 

10.3 

10.0 

-3.3% 

10 

9.8 

- 2 . 1 % 

> 200GeU 

9.2 

8.8 

-4.1% 

8.9 

8.7 

-2.5% 

> 300GeU 

7.6 

7.2 

-5.4% 

7.4 

7.1 

-4.3% 

lOOGeU < < 200GeU 

1.1 

1.1 

3% 

1.1 

1.1 

0 . 8 % 

200GeU < < 300GeU 

1.6 

1.6 

2 % 

1.5 

1.6 

6.3% 


- 17 - 























